Far-from-equilibrium state in a weakly dissipative model 
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We address, on the example of a simple solvable model, the issue of whether the stationary 
state of dissipative systems converges to an equilibrium state in the low dissipation limit. We 
study a driven dissipative Zero Range Process on a tree, in which particles are interpreted as finite 
amounts of energy exchanged between degrees of freedom. The tree structure mimicks the hierarchy 
of length scales; energy is injected at the top of the tree ('large scales'), transferred through the 
tree and dissipated mostly in the deepest branches of the tree ('small scales'). Varying a parameter 
characterizing the transfer dynamics, a transition is observed, in the low dissipation limit, between a 
quasi-equilibrated regime and a far-from-equilibrium one, where the dissipated flux does not vanish. 

PACS numbers: 05.40.-a, 02.50.Ey, 47.27.eb 
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One of the main challenges of non-equilibrium statis- 
tical physics is to understand which principles rule the 
description of nonequilibrium stationary states. Generic 
approaches, like linear response theory [l|, have been de- 
veloped for weakly driven systems (for instance gently 
stirred fluids, or systems subjected to small temperature 
gradients). In this situation, when the external forcing 
is very small, the system remains close to an equilibrium 
state, and the effect of the small drive is thus perturba- 
tive. A different situation, which also attracted a lot of 
interest, is that of dissipative systems, in which an energy 
flux takes place between scales rather than in real space. 
Energy is usually injected at large scales, and cascades 
down, through non-linear interactions, to smaller scales 
where it is dissipated. Examples include hydrodynamic 
turbulence 0, 3 1 , wave turbulence in fluids or plasma 
01 , ferrofluids 5 1 and vibrating plates d, 0| , fracture 
and friction 
foams 



12|. 



I , as well as granular materials [10|, |ll| and 



Dissipative effects are usually characterized 
by a dissipation parameter, like viscosity or inelasticity 
coefficients, which is zero in the conservative case. A nat- 
ural question is to know whether in the limit of small, but 
nonzero dissipation coefficients, the stationary state of 
the system becomes close to some equilibrium state to be 
determined. Qualitatively, one may expect that adding 
a tiny amount of injection and dissipation to a conser- 
vative system breaks energy conservation, but leads to 
small fluctuations around a given energy level selected 
by the injection and dissipation mechanisms. The sys- 
tem would thus merely behave as if it was at equilibrium 
at this energy, and a perturbative approach around an 
equilibrium state may be meaningful. 

Whether this scenario holds in general is certainly an 
open issue. Perturbative approaches around equilibrium 
states for dissipative systems have been proposed in the 
context of two-dimensional turbulence 13, 14, 15, 1^, 17, 
[3] , for which the flux of dissipated energy vanishes in the 
small viscosity limit [H. However, in other situations like 



three-dimensional turbulence [JJj or granular gases [l9| . 
the dissipated flux seems to remain finite for small vis- 
cosity, suggesting that the statistical state of the system 
does not converge to any equilibrium state. 

In order to give further insights into these issues, sim- 
ple solvable models may be helpful. Here, we study a 
stochastic transport model, namely a Zero Range Process 
2ol . 21 1, that describes in a schematic way the transfer 
of energy between scales, in the presence of injection and 
dissipation. To account for the specific organization of 
scale space, where small scales are much more "numer- 
ous" than large scales, we define our model on a tree ge- 
ometry. We show that depending on the energy transfer 
dynamics, the dissipative stationary state of the model 
converges in the weak dissipation limit either to an equi- 
librium state, or to a well-defined far-from-equilibrium 
state with a finite dissipated flux. 

The model is defined on a tree composed of M suc- 
cessive levels (see Fig. [1]); at any given level j < M, 
all sites have m > 1 forward branches that link them 
to level J ' + 1, so that the number of sites at level j is 
Sites are thus labeled by the level index j, and 



the index i — 1, . . . , within level j. The energy on 

each site (j, i) of the tree is assumed to take only discrete 
values proportional to an elementary amount Sq, namely 
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FIG. 1: Sketch of the model, iUustrating the tree geometry. 
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Sj^i = Tij^iEQ. Energy transfer proceeds as follows: an 
energy amount eg is moved, either forward or backward, 
along any branch between levels j and j ' + 1 with a rate 
Uj. In the absence of driving and dissipation, the dynam- 
ics satisfies detailed balance. Energy injection is imple- 
mented by connecting the site (1,1) to a thermostat of 



temperature Text = Poxt^ with a coupling frequency :/cxt- 
Dissipation proceeds through the random withdrawal of 
an amount of energy Sq at site (j, ?) with rate A^. The 
master equation governing the evolution of the probabil- 
ity distribution P{{nj i},t) reads: 



J 



dP 

dt 



+ ^ ^ ^'^J [Pi^j^i + 1> "j + l,(j-l)m+i - 1, {nq,r}, t) + P{nj^i ~ 1, + + 1, {n,,r}, t)] 

j=l i=i ;=i 

+ ^jP^^J^^ + 1' {'^9.'-}' ^) + ''<=xt [P{ni,i + 1, {n,,,}, t) + e-'3-'-'»P(ni,i - 1, {n,,,}, t)] . (1) 

j=i i=i 
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In this last equation, {uq j.} is a short-hand notation for 
all the variables that are not explicitely listed. Turning 
to the stationary state, we look for a steady-state prob- 
ability distribution of the form: 



Pst{{nj,i}) 



4nn 



(2) 



where (3j is an 'effective' inverse temperature (to be de- 
termined) associated to level j, and Z is a normalization 
factor. Inserting expression ([2]) of the stationary distribu- 
tion into the master equation ([T]) yields a set of equations 
to be satisfied by the parameters zj = exp(— /3j£o), for 
j = 2,...,M-l: 



j/j_i(zj_i - Zj) - mujizj - Zj+i) 



(3) 



with the boundary conditions 

i^cxt(e"'^°="^'' - zi) - mvi{zi - Z2) = Aizi (4) 

vm-i{zm-i - zm) = AmZm- (5) 

Note that these equations correspond to the local balance 
of the diffusive fluxes Vj(zj — Zj+i) and dissipative fluxes 
AjZj. In the absence of dissipation, namely if Aj = 
for all j, the equilibrium solution /3i = ... = Pm = 
/3oxt is recovered. To study the dissipative case, we need 



kj). The transfer rate Vj and the dissipation rate A-,- are 
balanced for a wavenumber kj = K given by 



K 



j/qX 1/(7-") 



(S) 



(7) 



which goes to infinity in the limit of small dissipation co- 
efficient D. Note that V[)/D is similar to the Reynolds 
number in hydrodynamics. For large K, we shall call the 
ranges kj <C K and kj ^ K the "inertial" and "dissipa- 
tive" ranges respectively. 

The solution of Eqs. Q, ^ and (O can be evaluated 
numerically. We are interested in the inertial range ba- 
havior, where energy transfer dominates over dissipative 
effects, so that we shall explore the solutions by varying a 
while keeping 7 fixed. We first compute the mean energy 
flux <I> injected by the reservoir. 



-/3oxteo 



(8) 



The flux $ is plotted as a function of a in Fig. [Ifa), for 
a broad range of values of the dissipation coefficient D. 
We observe a transition around the value a = — 1: for 
a < —1, $ — > when D ^ 0, while for a > —1, <& con- 
verges to a finite value in the small D limit. These two 
regimes are also clearly seen in Fig. [2] (b), (c), and (d) 



to choose a specific form of the frequency Vj and the by plotting the temperature Tj = p. 



dissipation rate Aj. 

v., = i^ofc? 



A generic form is the following: 



Dk]^ 



7>0, 



(6) 



where we have introduced a pseudo-wavenumber kj = 
m^~^, to map the tree structure onto physical space. Pa- 
rameters vq and D are respectively a frequency charac- 
terizing the large scale dynamics, and a dissipation coef- 
ficient. We impose the condition a < 7, so that dissipa- 
tion becomes the dominant effect at small scales (large 



= — £0 / In Zj as a 
function of In kj = (j — 1) In to. A first trivial observation 
is that K increases more rapidly when decreasing D for 
larger values of a. More interestingly, we observe that 
for a = —2 [Fig. Htb)], the temperature profile slowly 
converges to the equilibrium profile Tj'^ — (3~J^ when 
D ^ 0, while for a = 1 [Fig. [2Kc) and (d)], it converges 
to a well-defined nonequilibrium profile, which is linear 
for kj < K when plotting (3j as a function of In kj [see 
Fig- [2Ic)]. These results can be interpreted as follows. 
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FIG. 2: Numerical solution of Eqs. © and dS]) for 7 = 2 
and D — 10~^, with x given in (a); same symbols for all four 
figures, (a) Energy flux ^ as a function of a; tie full line is $0 
given in Eq. i)JOD . (h) Temperature Tj = versus Infcj for 
a = —2; the full line is Text, (c) f3j versus \nkj for a = 1; full 
line: defined in Eq. i|J4D . (d) Same data as (c) plotted 

as Tj, on the same window of Ink j as (b). Arrows in (b) and 
(c) indicate the value of K for each D. Other parameters: 
M = 50, m = 2, 1/0 = 1, i^cxt = 0.1, /3cxt = 1, £0 = 1- 



When the transfer mechanism is inefficient at small scales 
(a < —1), dissipative scales are not "feeded", so that en- 
ergy accumulates at large scales, generating an effective 
equilibrium. In the opposite case (a > —1), the transfer 
mechanism is efficient at small scales, thus "pumping" 
energy from large scales to dissipative ones. 

Most of the above behavior can be understood using 
a simpler form of the dissipation, which leads to analyt- 
ically tractable calculations. We assume that = for 
all j < M, leaving a nonzero dissipation rate Am only 
on the last level j = M of the tree. As a first step, we 
look for the solutions of Eqs. ([3]) and (|4]) with Aj = 0, 
j = 1,...,M — 1, without taking into account the dis- 
sipative boundary condition ([5]). We find a family of 
solutions, parameterized by the flux (f>: 



-/5cxt£0 



1 - $ 



1 



Bk' 



where j — 1, . . . , M, and with $0 and B given by 

Note that Zj smoothly converges to e~'^'=='*^° when $ ^ 0, 
and is a decreasing function of kj for <f> > 0. Interestingly, 



(9) 

(10) 
(11) 



Eq. ([9]) imposes an upper bound 'I'max on the flux $, 
which is determined by the condition zm > 0: 



1 
$0 



1 
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(12) 
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If a < —1, one finds for large M that <I>ma 
so that $max when M — s- cx). Accordingly, whatever 
the small scale boundary condition, the flux vanishes in 
the large size limit. In contrast, if a > —1, $max con- 
verges to $0 > in the large size limit Af — > 00; $0 goes 
to zero linearly with a when a — 1"*". 

We now use the dissipative boundary condition ([5]) to 
determine the precise value of the flux. Eq. (O states 
that the diffusive flux $ is equal to the dissipated flux on 
level j = M. This condition leads to $ = to*'^"-'^ Am^m, 
or using kM = rn^^~^, zm = *&/(^AfAM). Equating this 
value of Zm with that given in Eq. © for j = M, yields 
an equation for $, which is solved into 



1 



^mAa/ 



(13) 



Identifying kM with the value K defined in Eq. ([7]), we 
get Am = DK^ = vqK'^, and thus fci^AM = i/o^^+"- 
For a < —1, both /cm A a/ and $max are for large K 
proportional to so that their ratio is a constant. 

From Eq. (|13p , $ goes to zero as a finite fraction of the 
maximum flux $niax- Using K ^ 15-1/(7-")^ the flux <I> 
behaves in terms of the dissipation coefficient as $ I?^ 
when £) ^ 0, with /i = |1 + a|/(7 - a). From Eq. (O, 
Pj — > /3oxt when 13 ^ 0, as long as kj ^ K. Altogether, 
the effect of dissipation on the system may be considered 
as perturbative in the case a < — 1. The perturbation 
expansion is however singular, with a nontrivial exponent 
/i < 1. In the opposite case a > —1, ^max ^ $0 > 0, 
while AjmAa/ 00. Hence from Eq. $ is equal 

for large K (or small D) to the maximum flux ^max = 
$0, consistently with Fig. [^a) From Eq. ©, the 

temperature profile f3j = — e^^'^lnz^- converges to a well- 
defined nonequilibrium profile 



= 1(1 



lln% +/3cxt + IhiC, (14) 
£0 



with C = 1 + i'o{m — m^") / Vcxt- Although this profile 
has been obtained with a simplified version of the model, 
one sees on Fig. [Sfc) that (3j computed numerically in the 
original model also converges to (3^'^- Let us emphasize 
that PY"^ does not depend on parameters related to dissi- 



pation j24| , but only on parameters characterizing injec- 
tion and transfer. The temperature profile is continuous 
with a, namely Pj — > /3oxt when a — 1+. Note also 
that the coupling i^ext simply "renormalizcs" the inverse 
temperature /3ext; the low coupling limit corresponds to 
driving the system with a small effective temperature. 

To sum up, it turns out that an equilibrium approach 
to the stationary state of the present model in the weak 
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FIG. 3: (a) (Y) as a function of a, for D = 10"^, • Y is 
defined in Eq. i)J5|) . with here f{n) = n and g{k) — k^^ . 
The full line is the asymptotic value of (Y) for D ^ 0. (b) 
Normalized entropy S/Nk versus a, for the same values of D 
as (a). Other parameters: same as in Fig.\^ 



to N'^' if —1 < a < 0, while S is independent of Nk 
for a > 0. Hence the dissipative state characterized by 
Pncqiinj^i}) has a much lower entropy than the quasi- 
equilibrium state obtained for a < —1. In other words, 
the accessible volume in phase space is much smaller in 
far-from-equilibrium states than in equilibrium states. 

A major challenge for future work would be to char- 
acterize such asymptotic nonequilibrium states in more 
realistic models, and to understand their statistical fun- 
dations. The fluctuation properties, an d sp eciflcally the 
validity of Gallavotti-Cohen relations llOl. in these 
states, would be issues of great interest [26|. 



dissipation limit is meaningful if a < — 1. In this case, 
the probability distribution Pst({"-j,i}) converges, in a 
weak sense, to the equilibrium distribution of tempera- 
ture /3oxt- Although some deviations from equilibrium 
persist at small scale in the distribution for D ^ 0, the 
average value of "large scale" observables converges to 
the corresponding equilibrium value. Such "large scale" 
quantities, which are not sensitive to small scale details 
of the distribution, include observables Y defined as 

M / m^-^ \ 

y = J2Uk,) E /Koj ' (15) 

where /{n) is an arbitrary function, and g{k) satisfies 
kg{k) when fc — * oo. In the opposite case a > — 1, 
the probability distribution Pst{{njj}) converges when 
D — > to a well-defined nonequilibrium probability dis- 
tribution P„cq({nj,i}), given by Eqs. © and (HH) 0. 
The convergence of (Y) to {Y)cq for a < —1 and to 
(^)ncq for a > —1 is illustrated on an example in 
Fig. (Ha). Note that {Y)ncq depends on a, and thus on 
the energy transfer, while {Y)cq obviously does not. 

Other relevant statistical quantities are sensitive to the 
small scale details of the distribution, and have a more 
complex behavior. This is the case of the entropy 

^ = - E Pst{{n,,.})\nP,t{{n,,,}). (16) 

As for a < —1, Tj drops from Text to zero for kj fa K, 
it is natural to expect that the entropy S is proportional 
at large K to the number Nk of sites (j, i) with kj < K, 
namely Nk = mK/{m—l). This behavior is confirmed in 
Fig. [DJb) where the normalized entropy S/Nk is plotted 
as a function of a for different small values of D: S/Nk 
converges to a well-defined value for a < — 1, while no 
clear convergence is observed for a > —1. Such a re- 
sult can be understood in the framework of the simpli- 
fied model where Am = vqK" and Aj — for j < M. 
For large K, S is proportional to Nk if a < —1 and 
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